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Abstract 

Using the Riemann-Hilbert approach, we study the quasi-linear Stokes 
phenomenon for the second Painleve equation y^x ~ 2j/'' + xy — a. The 
precise description of the exponentially small jump in the dominant so- 
lution approaching a/x as |a;| oo is given. For the asymptotic power 
expansion of the dominant solution, the coefficient asymptotics is found. 

1 Introduction 

The classical Painleve equations are six particular ODEs yxx = R{x,y,yx) 
whose general solutions called the Painleve functions are free from movable 
critical points and can not be expressed in terms of elementary or classical 
special functions. The properties of the Painleve functions show the remarkable 
similarity to those of the classical transcendents of hypergeometric type like 
Airy, Bessel, Weber, Whittaker and Gauss functions, see monograph ||^ and 
surveys mentioned there. 

The problem of the asymptotic description of the Painleve functions near 
fixed critical points is one of the most interesting and important problems in 
the theory of the Painleve equations. The first result of such kind was obtained 
by Boutroux who studied the first and second Painleve equations, yxx = 
6j/^ -|- X and yxx = Sy"^ + xy — a, using direct methods based on the integral 
estimates. Boutroux has shown that, generically, the asymptotic as |a;| — > oo 
behavior of the first and second Painleve functions is given by the Weierstrafi 
elliptic function and by the elliptic sine, respectively. Periods of these elliptic 
asymptotic solutions depend on argx only. Along the rays argx = it -\- '^k for 
PI and arga; = ^fc for P2, k ^ Z, the general Painleve transcendent is described 
by the degenerate elliptic, i.e. trigonometric, functions. The asymptotic forms 
of the phase shift in the elliptic ansatz for Pf and P2 as argx belongs to the 
interior of the "elliptic" sectors are found in ^ , while the connection formulae 
for the phase shift for different sectors were obtained in ^ . Similar results for 
Painleve equations P4 and P3 were published in The nontrivial jump in 

the phase of the elliptic asymptotic ansatz across the rays of the trigonometric 
asymptotic behavior constitutes the nonlinear Stokes phenomenon. 
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In our terms, the quasi-linear Stokes phenomenon is the analog of the clas- 
sical, or linear, Stokes phenomenon. The latter was observed for the first time 
in the asymptotic analysis of the Airy equation, yxx = xy, and consists in the 
jump of the recessive, exponentially decreasing solution in the "shadow" of the 
dominant, exponentially large term. 

Boutroux has found that the Painleve equations admit 1-parameter so- 
lutions represented by the sum 

y ~ (power series) -I- (exponential terms) (1) 

of a power series with the leading term y^ix) which satisfies certain algebraic 
equation, e.g. 6j/q -|- x = for PI and 2j/q -|- xj/q = o- for P2, and of a transseries 
which is the sum of decreasing exponential terms. The latter can be obtained 
from the conventional perturbation analysis. Solutions of such kind are usually 
referred to as the truncated, separatrix, degenerate, quasi-stationary or instan- 
ton solutions. 

Since all equations of the perturbation analysis are linear homogeneous and 
non- homogeneous ODEs, there is no surprise to meet the quasi-linear Stokes 
phenomenon studying the degenerate solutions. 

The first, heuristic, description of the quasi-linear Stokes phenomenon was 
obtained in 0, There are known several attempts to find the analytic mean- 
ing of the quasi-linear Stokes phenomenon and to justify it using the resurgent 
analysis. The reader can find general settings of the Borel summation tech- 
nique for the nonlinear differential equations in the paper of Costin The 
technique is shown to be effective for certain nonlinear asymptotic problems. 
For instance, in Costin applied the resurgent approach to find the asymp- 
totic location of the first real pole of the degenerate solution of PI when the 
amplitude of its exponentially decreasing perturbation becomes large. 

Takei |Q proposed two ways of the use of the Borel summation for descrip- 
tion of the quasi- linear Stokes phenomenon in PI. His first idea is similar to one 
of jl^ and is based on the study of the Borel transform of the leading power 
series. Assuming that the Borel series has certain analytic structure on the 
boundary of its convergence, Takei obtained a simple formula which relates the 
jump of the recessive term in (Q) to the asymptotics of the Taylor coefficients 
of the Borel transform. However, the assumptions on the analytic structure of 
the Borel transform remain unjustified. In the second way, he studied the lin- 
ear differential equation associated to PI using the so-called "exact" WKB 
analysis. This potentially effective approach is based on the assumption of the 
Borel summability of the formal "instanton" series for the Painleve function 
which is still an open problem, see 

The coefficient asymptotics in the power series in (0) for the fifth Painleve 
equation P5 was found by Basor and Tracy JTzt up to a common constant factor. 
The approximate value for the latter was proposed in ijl^ using the numerical 
computations. The exact value of the Basor- Tracy constant as well as the proof 
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of their asymptotic formula were presented by Andreev and Kitaev in p8[ via a 
combination of the isomonodromy method and the Borel summation technique. 
We emphasize the importance of the paper as the first exact description 
of the relation between the coefficient asymptotics in the Painleve transcendent 
expansions and the quasi-linear Stokes phenomenon. 

In , Joshi and Kitaev studied the asymptotics of the leading power series 
coefhcients for the degenerate solution of PI and announced some results in the 
asymptotics of the transseries coefficients. The investigation is based on the 
use of the recursion relations and leaves a common factor for all the asymptotic 
coefhcients undetermined. The authors proposed the value of this factor noting 
that it could be found using the quasi-linear Stokes phenomenon and the Borel 
transform technique similar to presented in . 

The common idea of the approaches mentioned above is the use of the Borel 
transform, most frequently, to fix the solution having the formal power series as 
its asymptotics. The problem with this idea is that one needs to assume certain 
analytic properties of the Borel transform of the solution. 

Below, we study the quasi-linear Stokes phenomenon for the second Painleve 
equation P2, 

Vxx ^ + xy - a, a = const, (2) 

from the viewpoint of the isomonodromy deformation method jl^, |o[ |^. We 
pursue a two-fold goal, i.e. (a) to make an exact sense of the formal equation 
(0), and (b) to evaluate the asymptotics of the coefficients of the leading power 
series in (^. 

We stress that we do not use the Borel summation at any stage of our inves- 
tigation. Instead, we fix the solutions of P2 (j|) with the power series asymptotic 
behavior by the proper choice of the monodromy data of the associated linear 
equation. Unlike [T^ , our study is based on the direct asymptotic analysis of 
the Riemann-Hilbert problem via the nonlinear steepest descent method psf . 
In this approach, we associate the solution to the RH problem with a disjoint 
graph. This in turn allows us to separate effectively the exponentially small term 
from the power-like background, see Section 2. Using the outlined method, 
wc rigorously prove the existence of the solutions yi{x), 2/2(2:), yz(x) asymp- 
totic to a/x as 00 in the respective overlapping sectors argx e (— ^,7r), 
argx g (^,^), argx g (— tt, -1), and find the exponentially small differences 
Vk — yi, k,l = 1, 2, 3, which constitute the quasi-linear Stokes phenomenon. 

The collection of the functions yk{x), k — 1, 2, 3, forms a piece- wise holomor- 
phic function y(x) ~ ^ as — s- 00. The moments of this function immediately 
yield the asymptotics for the coefficients of the leading power series (Q) . 



2 Riemann-Hilbert problem for P2 

In Flaschka and Newell formulated the inverse problem method for P2. The 
method reduces the integration of (||) for |Rea| < \ to solution of the matrix 
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Riemann-Hilbert (RH) problem for an auxiliary function ^(A) which in its turn 
is equivalent to a system of singular integral equations. Using this method, they 
described the limit of small y and yx and obtained the determinant formulae 
for the rational and classical solutions of (|^). Later, the RH problem was ex- 
tensively studied by Fokas, Ablowitz, Deift, Zhou, Its, Kapaev, who considered 
the questions of existence and uniqueness of the RH problem solution ^ , 
and obtained various asymptotic solutions of the latter [ESl E4| . 



In accord with |20 , the Painleve function set is parameterized by two of 
the Stokes multipliers of the associated linear system below denoted by the 
symbols Sfc, k — 1,2,3. We begin with the case of decreasing as, x ^ +00 
Painleve transcendents corresponding to the special value S2 = 0. We observe 
that, for < argx < the RH problem graph can be transformed to the 
union of three disjoint branches. Using this fact, we define the dominant term 
yi (x) as the solution of the Painleve equation corresponding to the reduced RH 
problem with S2 = S3 = 0. The recessive, exponentially small term appears as 
the contribution of the disjoint branches in the RH problem graph. Thus the 
recessive term can be effectively distinguished from the dominant background. 

Below, we use the modification of the original RH problem of Flaschka and 
Newell [ pO[ which is more convenient for our purposes. Also, we restrict ourselves 
to the case a — ^ ^ Z. The excluded case of the half-integer a can be treated 
similarly, but requires a separate consideration due to logarithmic behavior of 
\E'(A) at the origin. 

Let us introduce the piece- wise oriented contour 7 — CUp^Up-Uf_-^^Jk which 
is the union of the rays 7^ = {A G C: |A| > r, arg A — ■|-|-^(A;— 1)}, k — 1,...,6, 
oriented to infinity, the clock-wise oriented circle C = {A € C: |A| = r}, and 
of two vertical radiuses p+ = {A G C: |A| < r, arg A = ^} and p_ = {A G 
C: |A| < r, argA — — §} oriented to the origin. The contour 7 divides the 
complex A-plane into 8 regions Qk, k G {^, r, 1, . . . , 6}. Qi and Qr are left and 
right halves of the interior of the circle C deprived the radiuses p+, p_. The 
regions Qk, k — 1, ... ,6, are the sectors between the rays 7^ and jk-i outside 
the circle. 

Let each of the regions ftk, fc = r, 6, 1, 2, be a domain for a holomorphic 2x2 
matrix function ^'^(A). Denote the collection of ^'fe(A) by 5* (A), 

*(^)lAen. = *'^(^)' = <^2*(A)a2. (3) 

Let 5*+ (A) and 'J- (A) be the limits of ^'(A) on 7 to the left and to the right, 
respectively. Let us also introduce the Pauli matrices 0-3 = ( J ) , o'2 = ( ^ ) ' 
(Ti = ( ° J ) and two matrices (t_|_=:("J),(t_ = (!|'"). The RH problem we talk 
about is the following one: 

i) Find a piece- wise holomorphic 2x2 matrix function 5* (A) such that 

^{\)e'^''' I, A ^00, e = i{^\^ + xX), (4) 
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and 

||*r(A)A""'"^ II < consi as A ^ 0; (5) 

ii) on the contour 7, the jump condition holds 

^+{X) = ^-{X)S{X), (6) 

where the piece-wise constant matrix 5(A) is given by equations: 
on the rays jk, 

S{X)\^^=Sk, S2k-i ^ I + S2k-io--, S2k = I + S2ka-+, (7) 

with the constants Sfe satisfying the constraints 

Sk+3 = -Sfc, si - S2 + S3 + S1S2S3 = -2sin7ra; (8) 

on the radiuses p±, 5(A) is specified by the equations 

Xep-: *((e^''*A) = *^(A)M, 

Xep+: -iiriX) = ^i{X)a2Ma2, (9) 

where 

M = ie*'^"'^Vi; (10) 

on the circle C, the piece-wise constant jump matrix 5(A) is defined by 
the equations 

*6(A) = *,(A)£;5-\ ■fi{X) = ^r{X)E, *2(A) - *.(A)£;5i, (11) 
^'3(A) = ^i{X)<j2E<j2S~\ *4(A) - ^i{X)<j2E<j2, ^'5(A) = ■i'i{X)a2Ea2S4, 

where the unimodular constant matrix E satisfies the equation 

ES1S2S3 = a2M-^Ea2. (12) 

Because the asymptotics of 5* (A) as A ^ 00 is given by 

F(A) := vI/(A)e''- = I +l{-'-Ra, + ^a,) +0{^), D^yl-y^~xy^ + 2ay, 

(13) 

the solution y{x) of the Painleve equation can be found from the "residue" of 
Y{X) at infinity, 

y = 2 lim AYi2(A) = 2 lim Ar2i(A). (14) 

A — 'oo A — *oo 

The equation (^4|) specifies the Painleve transcendent as the function y = 
f{x,a,{sk}) of the independent variable x, of the parameter a in the equa- 
tion and of the Stokes multipliers Sk (in the generic case, the connection matrix 
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E can be expressed via Sk using the equation (|T^) modulo the left diagonal mul- 
tiplier). Using the solution y — f{x,a,{sk}) and the symmetries of the Stokes 
multipliers described in we obtain further solutions of P2: 

y = -f{x,-a,{~Sk}), y ^ f{x,a,{s^}), 

y = e*¥"/(e^¥"x,a,{sfe+2„}), (15) 

where the bar means the complex conjugation. 
2.1 Asymptotic solution for S2 = 

Let us consider the RH problem above where S2 = assuming that ^ cx3 
in the sector argx e [0, Equations (||) imply that S5 = as well, and the 
function Y{X) defined in ( p^ ) has no jump across two rays arg A = ^, |A| > r. 
The remaining parameters si and S3 satisfy the constraint 

si + S3 = —2 sinTra, 

and the connection matrix E due to (^2|) has the form 

Our first step in the RH problem analysis is elementary and consists of 
the formulation of the equivalent RH problem for the piece-wise holomorphic 
function Y{X) = 'J'(A)e''(^)'"3 : 

i) r(A)^/ as A ^00, \\Y{\)\-°"'''\\ < const as A ^ 0; 
a) r+(A) =y_(A)G(A), G(A) = e-^'"^5(A)e^'"^ A e 7. (18) 

If S{X) = I + sa± then G{X) = I + se^'^^a±. On the second step, we 
transform the jump contour 7 to the contour of the steepest descent for the 
matrix G(A) — /. 

Let us denote by 7+ the level line lm6(X) = const passing through the 
stationary phase point A+ = ^x^^^ and asymptotic to the rays arg A = 
This is the steepest descent line for e^^. Similarly, let us denote by 7_ the level 
line passing through the stationary phase point A_ = —^x^^"^ and asymptotic 
to the rays arg A = — ^, This is the steepest descent line for e~^^. 

For arg a; G (0, the level lines i+ and £- emanating from the origin and 
asymptotic to the respective rays arg A — ^ and arg A = ^ are the steepest 
descent lines for e^^ and e^^^, correspondingly. 

For argx = 0, the level lines £+ and ^_ are the segments connecting the origin 
A = with the respective stationary phase points A+ and A_. For x > 4r^, these 
lines contain the radiuses /?+ and p_ . 

It is convenient to consider the following equivalent RH problems for ^'(A): 
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(16) 
(17) 



for arga; S (0, the jump contour is the unfon of the level lines 7+, 7_, , I- 
oriented from the left to the right, and of the circle C = Cr U C; divided by i± 
into the right Cr and left Ci arcs both clockwise oriented. The jump matrices 
are as follows: 

Ae7+: ^(A)-5^^= (^_^^^ J), Ae7-: ^(A)-56- (J 7' 
A.^,, |A|>.: 5(A) = 5,^3= ;)^5„ 
A e |A| > r: 5(A) = Sg-i^ri = J -2sin™ j ^ 

A e |A| < r: 5(A) = G-2.M~^a-2., A S |A| < r: 5(A) = M, 

A G Cr- 5(A) = i;, AG Q: 5(A) = G2Ea2- (19) 

The jump contour for the RH problem ( |l9|) is decomposed into the disjoint 
sum of the lines 7+, 7_ and 70 = £+U£_UC For the boundary value arga: = 0, 
the lines ^+ and i- emanating from the origin pass through the stationary phase 
points and partially merge with the lines 7+ and 7_. Thus the corresponding 
RH problem graph can not be decomposed into a disjoint union of the steepest 
descent contours except for a G Z when the jump across 70 can be eliminated. 
The particular case a = was studied in |^. The general case is described 
in the following way: 

for argx = 0, the jump contour is the union of the level lines 7+, 7_ oriented 
from the left to the right, of the level lines ^+ and £_ oriented from A_ to A_|_, 
and of the clockwise oriented circle C — Cr^Ci. The jump matrices are as 
follows: 

S^^, ReA < 0, ^ ^ „, . Qi^^ _ / S^^, ReA < 0, 



'^'^-<SVCX>C ^^^-■■'^') ,5e,ReA>0, 



A G £+, |A| > r: 5(A) = 5i53 



1 
— 2sin7ra 1 



A G |A| > r: 5(A) = S^^S^^ = ' ^ 2sm™ 



1 

A G £+, |A| < r: 5(A) = craAf-Va, A G |A| < r: 5(A) = 

A G Cr- 5(A) =E, AG Ci: 5(A) = a2Ea2. (20) 

The problem ( po|) is the limiting case argx — * of (^9|). 
As arga; G [0, -1], introduce the reduced RH problem {s2 = S3 = 0) for the 
piece- wise holomorphic function $(A) on 70 = ^+ U U C oriented as above: 

i) $(A)e^'"^^/, A ^00, ||$(A)A""''^|| < coTisi, A ^ 0, 
ii) $+(A) = a>_(A)5(A), AG 70, (21) 

where the jump matrix 5(A) is described in ([l9|). 
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Theorem 2.1 If a - \ ^ Z and argx € [0, f-] while \x\ is large enough, then 
there exists a unique solution of the RH problem j^l 



Proof. Since det S{X) = 1, we have det = det and hence det $(A) is 
an entire function. Furthermore, because of the normahzation of <I>(A) at infinity, 
det«'(A) =. Let $ and $ be two solutions of (|2l|). Then x(A) = l>(A)$"^(A) is 
a rational function of A with the only possible pole at A = 0. However, due to 
the boundedness of <^{X)X'°"'^ and $(A)A-"'^=' as A ^ 0, we have x(A) bounded 
at A = and, using the Liouville theorem, x(A) = const. Thus, because of the 
normalization of $ and $ at infinity, <I>(A) = ^(A), and the solution is unique. 

To prove the existence of the solution <i>(A), introduce an auxiliary hmction 



1 IW 1 

-1 1 / I 1 



(22) 

where the prime means differentiation w.r.t. z and, for a — i ^ Z, 



4feA:!r(a + i + /c) ^ 2^ ""^^ ^' 

°° PI'S _ n-'lp-l^^+^fe 

Here Jt/(z) denote the classical Bessel function |^5|. It is worth to note that the 
function <t'^{z) satisfies the linear differential equation 

<|0^(^3_^^2)$°. (24) 

As it is easily seen, if \z\ < const, then 

\\^°{z)z-°"'''\\ < const. (25) 
Using the properties of the Bessel functions, we find that the products 

$i(z) = $*'(z)^, ^2{z) = ^i{z)Si, $3(^)^*2(^)^2, (26) 



where 



A_/l2sin7ra\ Q _ f 1 
^i^lo 1 j' *2-U2sin7ra 1 



satisfy the asymptotic relation 



<^'fe(2) = (/-gfTl+0(l))e^-^ z^oo, argzG (7r(fc~|),^(fc+^)). (27) 
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We observe also the symmetries, 

a2$°(e"z)f72 = $"(z)M, Af = ie"""Vi, (T2$fc+i(e*"z)a2 = $fc(z), (28) 
and the relation 

(29) 



D 




Introduce the graph 70 = R U C consisting of the real axis oriented from 
the left to the right and of the clockwise oriented circle C — {z ^ C: \z\ — f} 
divided by the real axis in the lower Cd and upper Cu arcs. This graph divides 
the complex z-plane into four regions: which is the exterior of the circle in 
the lower half of the complex z-plane, Ct^ which is the exterior of the circle in 
the upper half of the complex z-plane, fid which is the lower half of the interior 
of the circle and Cl^ which is the upper half of the interior of the circle. Define 
a piece- wise holomorphic function $(z), 

$3(2), z e Cls, /3Q>, 
$°(z)i?, zeUd, ^"""^ 

CT2$°(e-'"z)i^(T2, Z £ 

By construction, this function solves the following RH problem: 

i) l>(z)e"^'"^ ^ /, z^oo, II $(z)z""'"3 II < co?isi, z ^ 0; (31) 
a) on the contour 70, the jump condition $+(z) = $_(z)5(z) holds true: 

z>f: S{z)^S2=(^_^^^^^ 

z<-f: ^(z) = 5-=(j -27™)=5_, 

-f<z<0: S{z)^M, 0<z<f: S'(z) = CT2Af"V2, 

zeCd-. S{z)^E, ze(7„: S{z) ^ a2E<j2. (32) 

Therefore the function $(z) has precisely the jump properties of the function 
$(A). To find $(A) with the proper asymptotic condition at infinity, let us 
replace the jump contour 70 by the contour obtained from 70 using the mapping 

z{\)^-e{X)^~ixX~i-X^. (33) 

3 

For |A| < R < ||x|^/^, the mapping ( ^3| ) gives us the holomorphic change of the 
independent variable. Introduce the piece- wise holomorphic function $(A), 

1,,^, f$(z(A)), |A|<i?, .34. 
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We look for the solution of the RH problem ( pT| ) in the form of the product 

$(A) = x(A)l>(A). (35) 

With this purpose, consider the RH problem for the correction function xW 
on the union £ of the clockwise oriented circle C of the radius R divided by the 
level lines £+ and £_ into the left £i and right £r arcs and of the outer parts of 
the lines £+, £_: 

i) XW I, A ^ oo; 

ii) x+(A) -x-(A)i/(A), \e£, (36) 



A e £+, |A| > R 
A G |A| > R 
|A| = i?, A e C 
|A| AeA 



i^(A) = e-^'"^S'+e^'^^ 

i/(A) = <|.2(z(A))e«-^ 
i/(A) - <i>3(2(A))e«--\ (37) 

Taking into account the equations (|3|), 5+ = 5*2, S- — S^^ , it is easy to 
see the continuity of the RH problem at the node points ^+ n £ and £- n £. 
Furthermore, the jump matrix i?(A) satisfies the estimates 

where the concrete value of the constants Cj, j = 1,2, is not important for us. 
Because we may take R = c|a;|^/^ with some positive constant c, the solvability 
of the RH problem (^) and therefore of (^l|) is straightforward. Indeed, consider 
the equivalent system of the non-homogeneous singular integral equations for 
the limiting value X~(A), i.e. 

or, in the symbolic form, — I ~^Kx~ ■ Here A_ means the right limit of A on 
£, and K is the composition of the operator of the right multiplication in H — I 
and of the Cauchy operator C_. An equivalent singular integral equation for 
ip~ X^ ^ I differs from ( |39| ) in the nonhomogeneous term only, 

il:-=KI + K^-. (40) 

Consider the integral equation ( |40| ) in the space H^{£) of the functions -0" 
such that and its distributional derivative both belong to L'^{£)- Since 
H — I is small in H^{£) for large enough \x\, and C_ is bounded in H^{£) pof , 
then ||_fi'||/fi(^) < c|a;|~^/^ with some constant c, while / — 7^ is invertible in 
H^{£). Because KI e H^{£), equation (EH) for is solvable in H'^{£), and 
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the solution xW of the RH problem ( ^ ) is determined by ^ (A) using the 
equation x = I + KI + Kij;^ . □ 

In accord with the said above, the function is given by the converging 
iterative series, x~ = X]^o Thus, for large enough \x\ and argx € [0, 

the Painleve function yi{x) corresponding to 82=3^ — 0, si = — 2sin7ra, is 
given by (|l4|), 

2/1(2:) = 2( lim x{X)X)^^ = -^^ /(x-(C)(i?(C)-/))i2^C- (41) 

A — >oo lib J 1^ 

Since x"(C) = + C(C"^) and \C,\>R^ c\x\^''^ , we find 

yi{x) - / (i?(A) - /)i2 d\-{I + 0{x-^/^)). 



Using here the expressions for the jump matrix H(\) from (|36D, the asymptotics 
(p7|), and the definition (^), we find the leading asymptotic term of yi(x), 

2/i(x) = ^ + 0(x-3/2), arga;e[0,^]. (42) 



Let us go to the case of the nontrivial S3 described by the RH problem (|19[). 
We look for the solution ^'(A) in the form of the product 

*(A) =X(A)$(A), (43) 

where ^(A) is the solution of the reduced RH problem (^T|). Thus we arrive at 
the RH problem for the correction fmiction X(A), 

i) X{X)^I, X^oo, 

a) X+{X) = X^iX)n{X), Ae7+U7_, (44) 
A £7+: n{X}^^^{X)S^^^z\\), 
Ae7_: H(A) = $_(A)56$I^(A), 

where 3>_(A) denotes the right limit of $(A) on 7+ U 7_. 
The estimate for the jump matrix on 7±, 

||W(A)-/|| < ce^m'^' -^s^)e-^\'=\''"\^-^"\\ (45) 

where c is some constant and Aq = ±|a;^/^ is the stationary phase point, yields 
the estimate for the norm of the singular integral operator K. in the equivalent 
system of singular integral equations, X_ = / + /CX_ , 

\\IC\\m^+u^^) < ce-il-l'''-^(i-s-). (46) 

If \x\ is large enough and argx E [0, ^ — e], e > 0, then the operator K. is 
contracting and the system X_ = I + JCX_ is solvable by iterations in {j^ U 
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7-), i.e. X- — Yl°^=a ^"^X^ = I + K,I + 0{e^^^ ' ). However, to incorporate 
the oscillating direction argx = -I in the general scheme, we use some more 
refined procedure. 

Theorem 2.2 If S2 — Q and \x\ — s- oo, argx € [0,^], then the asymptotics of 
the second Painleve transcendent is given by 

y = y,{x, a)~^ x-^l'e'i-"\l + ^ o{s\e-^^-''\ (47) 

where yi{x,a) ^ ^ is the solution of the Painleve equation for S2 — S3 = 0, 
si = —2 siuTra. 

Proof. Using the asymptotics of $(A) at infinity (pT|), we find the asymptotics 
of the jump matrix 7i(A), 



^' ^ I- S3e-^'^cr++0{ij-e-^^), A G 7_, ^ ' 



where '0~ = {I — K)~^KI introduced above is holomorphic in a neighborhood 
of 7± if argx £ (0, |] and ip~ G H^{l+ U 7_) if argx = and satisfies the 
estimate 

\\rm<c{\x\''^ + \\\r\ x&i±, 

with some constant c which value is not important for us. 
Consider the triangular matrices P(A), Q(A), 

S3 f ^26(0 _ dC 



*3 f ,-26(0. '^^ 



which solve the triangular RH problems: 

i) P(A) ^ / as A ^ cx), and P+(A) = P-(A)(/ - sse^^cr.) as A e 7+, 

ii) (5(A) -> / as A -> 00, and Q+{X) = Q-{X)iI ~ S3e~'^^a+) as A e 7_. 
We will look for solution X{X) of the RH problem (^) in the form of the 

product 

X{X) = Z{X)Q{X)P{X). (50) 
The correction function Z{X) satisfies the RH problem 

i) Z{X) I, X 00, 

It) Z+{X) = Z_{X)V{X), Ae7+U7-, (51) 

where 

V{X) = Q(A)P_(A)7i(A)(/ + S3e^'a^)Pl\X)Q-\X), X E 7+, 
ViX) = Q_(A)P(A)K(A)P-1(A)(/ + s3e-^'a+)Qz\X), A e 7_. (52) 
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Using the first of the relations in (48), we find the estimate for the jump matrix 
V{X) on 7+, 

1/(A) =/ + O(V'-e2«)+O(0(A)e2«), A e 7+, (53) 

where <^(A) = e"^^*^''^ is square integrable and holomorphic in a neigh- 
borhood of 7+. This function satisfies the estimate 

|0(A)| < c(|x|i/2 + |A|)"ie-3l-l'^'^-(i"g-), A e 7+, 

with some constant c which value is not important for us. For the product 
PJiP^^ on 7_, we compute the expression 

P{X)n{X)P-\X) =1- S3e-^'a+ + 0((^(A)e-28) + 0(^-e-2«), A G 7_, 

where (p{X) = J^^ s'^^^''^-^^ is square integrable and holomorphic in a neighbor- 
hood of 7_. This function satisfies the estimate 

\ip{X)\ < c(|x|i/2 + |A|)-ie-il"l'^'^°''(i A e 7_, 
with some constant c. Thus, 

F(A) =/ + C'(V"e-2S)+C'(v3e-^''), A G 7_. (54) 



Our next steps are similar to presented in the proof of Theorem 2.1. Consider 
the system of the singular integral equations for Z-{X) equivalent to the RH 
problem (|l]), = I + KZ_. Here the singular integral operator K is the 
superposition of the multiplication operator inV — I and of the Cauchy operator 
C_. Because the Cauchy operator is bounded in -^^(7+ U 7_), the singular 
integral operator K for large enough |a;|, argx € [0, f], satisfies the estimate 

||K|U^(,,U7-) <c|:rrl/2e-f|.|^^^cos(|arg.)_ (55) 

Thus equation C_ = K/ + K(^_ for the difference C_ := Z_ — / is solvable by 
iterations in the space ^^(7+ U7_). Solution of the RH problem ds^) is given by 



the integral Z ^ / + K/-KKC-. Thus, using (14), (M), IM) and the asymptotics 



of P, Q from {\m, we arrive at (47). □ 



2.2 The decreasing degenerate Painleve functions 

Applying the second of the symmetries ( p^ ) to (^), we obtain 

Theorem 2.3 If S2 =0 and \x\ 00, argx G [—-1,0], then the asymptotics of 
the second Painleve transcendent is given by 

y = y,ix, a) + ^ x-^l'e~^^^"\\ -I- 0{x-^l'')) + 0{s\e-'^^''^\ (56) 



•where 2/3(0;, a) = y\(x, a) ^ ^ is the solution of the Painleve equation for S2 
Si = 0, S3 = — 2sin7ra. 
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The solutions yi{x, a) and 2/3 (x, a) = yi(x, a) are meromorphic functions of 



j; G C and thus can be continued beyond the sectors indicated in Theorems 2.2 



and 2.3. To find the asymptotics of the solution 2/3(0:, a) in the interior of the 
sector arga: € [0, |], we apply (p7[). Similarly, we find the asymptotics of the 
solution 2/1 (a:, a) in the interior of the sector arga; g [~f jO] using ([5^). Both 
the expressions imply that, if |a;| — > 00, argx G [^f i fli 

/ \ I \ .siuTra _! 3/2 -i/2\\ , -ix3/2 

ya a;," - 2/1(2^:" = * — '/*e (1 + 0(a; + e 3^ . 57) 

Remark 1. The same asymptotic relation can be obtained using the analysis 
of the RH problem for the ratio ^'(A)^'~^(A) where ^(A) is the solution of the 
RH problem for S2 = si = while ^'(A) is the solution for S2 = — 0. 

Besides ?/i(a;, a) and 2/3(0;, a), let us introduce the solution 2/2(2;, a) corre- 
sponding to the Stokes multipliers si = S3 = 0, S2 = 2sin7rQ:. Due to the 
second of the symmetries (|l5|), 2/2 (x, 01) = y2{x, a). 

The last of the symmetries (|5|) yields the relations 

2/1 (x, a) = e*^2/2(e'^a;,a) = e~*^2/3(e"*^x, a), 
2/2(2;, a) = e*^2/3(e*^a;,a) = e"*^2/i(e"*^x, a), 
2/3(x, a) = e*^2/i(e'^2;, a) = e"*^2/2(e"*^x, a), (58) 

which imply the decreasing asymptotics of the functions ykix, a) in the following 
sectors: 

2/1 (x, a) , \x\^oo, argxG[0, — , 

X 6 

I . OL 2ti Att 

2/2(x,a) , |x| ^ 00, argxe[— ,— J, 

X 6 6 

(X 27r 
2/3(x,a) -, |x| -> 00, argxe[ — ^,0]. (59) 

The same symmetry (^5|) applied to the equation (^7|) yields the differences 

arga; e [-7r,--J: 

r \ / \ siuTra _-|/4 2 3/2 ^-r^/^^^ 

TT 

arga; ^ [g^^^]: (60) 

/ \ , \ siuTra _ma 2^3/2 _i/9\n It^/^, 

2/i(x,a) -2/2(x,a) = — -^x ^'^e^'' (1 + ©(a; ^Z^)) + ©(ea^ ). 



Equations (|5^) , (|60|) constitute the quasi- linear Stokes phenomenon for the sec- 
ond Painleve equation. 
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Finally, applying the last of the symmetries ( p^ ) to (p7|), (|5^), and using 
8|) we find the asymptotics of the degenerate Painleve functions: 

i) if si = 0, then S3 — S2 = —2 siuTra, and 

argx e hvr, -|^]: = y,{x, a) ~ ^ x-^l^el^''\\ + 0{x-^l^)), 
argxeh^,-^: = y3(x, «) - x-^^ef-^^^ (1 + 0(:,-V2)). (g^) 

ii) if S2 — 0, then si + S3 = —2 sinTra and 

argx e [-^,0]: y = 2/3(2;, a) + x-V4e-|-^'=(l + 

argx e [0, J: y = a) - x~^l \-"^-''\\ + ©(^-i/^)); (62) 

iii) if S3 = 0, then si — S2 = —2 sinvra and 

arg:r G [^,^]: y = 2/2(:.,a) - x-^l\l^'"\\ + ©(x-^/^)). (63) 

3 The power expansion of the degenerate solu- 
tions and the coefficient asymptotics 

Using the steepest descent approach, cf . Q , we can show the existence of the 
asymptotic expansion of yk{x,a), k ~ 1,2,3, in the negative degrees of x^^^. 
Further elementary investigation of the recursion relation for the coefficients 
of the series allows us to claim that the asymptotic expansion for any of the 
decreasing degenerate solutions has the following form: 

00 

yo{x, a) = - V + 0{x~n, (64) 

X ^-^ 

n=0 

where coefficients a„ are determined uniquely by the recurrence relation 

n 

ao = 1, a„+i = (3n + l)(3n + 2)a„ - 2a^ ^ akaiam- (65) 

kJ,m = Q 
k + l + m = n 

The initial terms of the expansion are given by 

af^ , 2(1 -g^) , 4(10-13«^ + 3«'t) , 
yo{x,a)^-[l + + + (66) 
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— ^ (280 - 397a^ + 129a'* - 12^^) 
16 



r (15400 - 22736a2 + 8427a'' - 1146aS + 55a*) +0(a;"i^)|. 

For a = 0, the recurrence (^5|) is exactly solvable, 

r(3n) 



3"-ir(n)' 



a = 0. (67) 



Our next goal is to determine the asyniptotics of the coefficients a„ in (^) as 
n ~* oo for arbitrary a. With this purpose, let us construct a sectorial analytic 
function y{x), 

(y3{x,a), arga; e (-^,0), 
y{x) ^ { yi{x,a), arga;G(0,%), (68) 
[y2ix,a), aigx£{^,^). 

The function y{x) has a finite number of simple poles. Therefore y{x) is bounded 
for |a;| > p and has the uniform expansion (64) near infinity. Let y^^\x) be a 
partial sum 

N-l 

yW(x,a) = - Va„x-3", (69) 

n=0 



and v^^''{x) be a product 

oo 

= x'^'iyix) - yW(x, a)) = - ^ a„+jvx-'" + 0{x-°°). (70) 



n=0 



Because x^^j/'-^-' (x, a) is polynomial, the integral of ?;'^^^(x) along the circle of 
the radius |a;| = p satisfies the estimate 



v^^Hx)dx <p'^^(f |y(a;)|d/ < 2V^+'max|y(x)| =Cp3^+* 

x\=p ■J\x\=P 



(71) 

On the other hand, inflating the sectorial arcs of the circle = p, we find that 

v^^\x)dx=i> v^^\x)dx+ 

+ r x^''{y^-y^)dx+ r x'^iy2^yi)dx+ x^"" {y^-y^) dx. (72) 

J p Je' 3 p Je ' 3 p 

Because v'^^^x) — ^o-n + 0{x^'^), the first of the integrals in the r.h.s. of ( [72[ ) 
is computed as follows: 

v'-'^^{x)dx = 2TTiaaN + 0{R-^). (73) 

\x\=R 
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Last three integrals in (|7^) are computed using (pT\), (|60|J 

x^^{y,-y3)dx+ x'^iy2-yi)dx+ r x^"" [ys ~ dx ^ 

= 3 rx3^(2/i-2/3)d2: = -3z^^ r x^''~-^e-i-"'\l + 0{x-''^))dx^ 
Jp V"" Jp 

^^g^sin™^3^2^_ir ^2Ar+l)(i + o(Ar-V 

Thus, letting R = cx), we find that the asymptotics as iV oo of the coefficient 
OAT in (|6j) is given by 

a^ = ^^(3/2)2^+*r(2iV+i)(l + 0(iV-i/3))+0(C^), N oo. (74) 

For a = 0, this equation is consistent with equation (|67|). For a G Z\{0}, 
the asymptotics ( |7^ ) is reduced to aat = 0{C^) which is consistent with the 
expansion of a rational solution of P2 determined by the triviality condition 
Si = S2 = S3 = and a G Z. 
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